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Abstract. We study the evolution of non-linear circularly polarised Alfven waves by solving numerically the 
time- dependent equations of magnetohydrodynamics (MHD) in one dimension. We examine the behaviour of the 
waves and find that different physical mechanisms are relevant in different ranges of /3. In a low j3 plasma the wave 
may undergo a parametric decay. This is because the wave excites a density enhancement that travels slower than 
the wave itself and thus interacts with the wave. When /3 > 1 the density enhancement does not interact with 
the wave and no decay takes place, instead the Alfven wave is reflected against the density enhancement. The 
reflection zone propagates with the speed — va- Because of that the magnetic flux is conserved which results in an 
amplification of the oscillating magnetic field by a factor i . We find that n depends on (3, and that in particular 
it is < 1 for values of ft ~ 1 and > 1 for /3 > 1. We discuss the relevance of these mechanisms to the acceleration 
of the solar wind, and the triggering of MHD turbulence in the polar wind region. In particular these simulations 
can explain the presence of inward propagating Alfven waves in the solar corona. 

Key words. MHD - waves - solar wind - stars : mass-loss. 



1. Introduction 

The expanding coronal model proposed by Parker (1958) 
predicted the existence of the solar wind. Observations 
by space probes showed later that there are two forms of 
solar wind: the slow solar wind (up to 400 km/s) which 
fits the Parker model, and the fast solar wind (up to 800 
km/s) which emanates from the coronal holes, regions in 
the solar corona with open magnetic field lines and low 
density. 

Large amplitude, low frequency, Alfven waves have 
been observed in the solar corona for over 30 years (e.g. 
Belcher & Davis 1971). During the last decade Ulysses has 
provided plasma and magnetic field measurements that 
have allowed extensive investigations of the behaviour of 
Alfvenic turbulence in the high-latitude solar wind. The 
data shows a strong correlation between the fluctuations 
in velocity and magnetic fields (Smith et al. I1995|l , and re- 
vealed the presence of both inward and outward directed 
Alfven waves (e.g Bavassano et al. 120000 . While the out- 
ward going waves are expected to be generated at the 
coronal base, the source of the inward going waves in the 
solar wind is not yet understood. 

Alfven waves play a crucial role in several models for 
the acceleration of the fast solar wind (e.g. Leer et al. 
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1982J. An attractive feature of the Alfven waves is that 
they can propagate over vast distances since they are in- 
compressible to lowest order, and therefore do not dissi- 
pate easily. However, dissipative damping is required at 
some point to avoid too high wind velocities (e.g. Holzer 
et al. ESSJ- 
In two dimensions phase mixing due to a transverse 
gradient in the phase velocity (e.g. Heyvaerts & Priest 
1983} can lead to a strong damping of the wave. This 
mechanism is not available in one dimension though, and 
one has then rather to consider the nonlinear coupling of 
the Alfven wave with other modes (e.g. Wentzel I1973fl . 
A close examination of the properties of Alfven waves 
shows that a linearly polarised Alfven wave is compress- 
ible to second order though, because the magnetic pres- 
sure, B\j2[iQ, is modulated on half the wave length of 
the Alfven wave itself (Alfven & F alt hammar 1 1 9 6 5|l . This 
opens up the possibility that a high-amplitude Alfven 
wave can steepen, which was demonstrated analytically 
by Cohen & Kulsrud l|1974f) . As the wave steepens, it 
forms current sheets at the nodes of the fluctuating mag- 
netic field. This effect has been studied in numerical sim- 
ulations by Boynton & Torkelsson l|1996[l and Ofman & 
Davila (JMZ}. 

In a circularly polarised Alfven wave, on the other 
hand, the magnetic pressure is constant along the wave, 
which is the physical reason why it is an exact solution 
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Fig. 1. Low amplitude magnetohydrodynamic waves 
propagating through a homogeneous medium of low (3 
(Run la), a Ap/p = ^zBi versus z/X at t/P = 4.2. b 
B x /B z (solid line) and v x /va (dotted line) versus z/X at 
the same time. 



of the nonlinear MHD equations. However, the wave can 
still decay via a parametric instability (e.g. Sagdeev & 
Galeev E!69i), which is usually less important than the 
wave steepening for a linearly polarised Alfven wave. In 
this instability a forward propagating Alfven wave in the 
presence of a density fluctuation generates a forward prop- 
agating acoustic wave and a backward propagating Alfven 
wave. 

Parametric decay has many applications in plasma 
physics and astrophysics. In coronal physics it has been 
proposed to be a possible mechanism to trigger MHD tur- 
bulence in regions with relatively smooth density profiles 
like the polar wind region and to account for the small 
compressible fluctuations in the solar wind (Goldstein 
I1978|l . It is also a candidate for generating the inward 
propagating Alfven waves in the solar wind (e.g. Tu & 
MarschEHEni) 

The parametric decay of circularly polarised Alfven 
waves has been studied both analytically (e.g. CohenEZHl 
Derby IT373I Goldstein HTMl Jayanti & Hollweg IT553jl 
and numerically by several groups (e.g. Del Zanna et al. 
ISnnil Malar a et a l. 127171711 Pru neti & Velli 115571 Ghosh 
& Goldstein IT5541 Ghosh et al. HTMl Umeki & Terasawa 
Except for the work by Pruneti & Velli l|T557)> these 
studies have been restricted to Alfven waves in homoge- 
neous media. Most of these numerical simulations have 
studied how Alfven waves of different frequencies interact 
and generate turbulence, and determine the properties of 
fully developed turbulence. For these purposes it is appro- 
priate to use a model with periodic boundary conditions 



and introduce Alfven waves that extend over the entire 
grid through the initial state (Malara et al. 120001 and Del 
Zanna et al.EEDJ- 

On the contrary we will study how a fresh Alfven wave 
falls prey to the parametric instability, which gradually 
converts it into a backward propagating Alfven wave and 
an acoustic wave. To do this we drive the Alfven wave 
on one of the boundaries of our model. Since there is no 
Alfven wave in the interior initially we are able to study 
how processes at the propagating wave front eventually 
lead to the breakdown of the entire Alfven wave further 
upstream, a phenomenon that cannot be found using peri- 
odic boundary conditions. Our approach is similar to that 
of Boynton & Torkelsson l|199fi|) . who studied linearly po- 
larised Alfven waves. Since the processes that we study 
occur gradually it is important to use a grid that extends 
over many wave lengths to capture them, which limits us 
to using a one-dimensional model. The need for an very 
extended grid becomes clear if one compares the works by 
Torkelsson & Boynton iflTMjl and Ofman & Davila l|T557|) . 
Torkelsson & Boynton found a much stronger damping in 
their one-dimensional model than Ofman & Davila did in 
their two-dimensional model, but most of the damping oc- 
curred beyond the outer boundary of the Ofman & Davila 
model. 

A highly simplified numerical model such as the one 
that we present in this paper cannot provide a realistic 
representation of the physics of Alfven waves in the solar 
wind, however it can still be useful since it allows us to 
study a limited number of physical processes in detail. 
The understanding that we gain from this can then serve 
as a guide in interpreting some aspects of more complex 
numerical simulations, such as the ones carried out by 
Tsiklauri et al. (T2"717l2"|) and Tsiklau ri & Nakariakov <12HTJ2")| 
and Laveder et al. I|2()02al l2002bf) . 

The plan of the paper is the following: In Sec. 2 we 
describe the basic MHD equations and review the proper- 
ties of the parametric decay. Our models and results are 
presented in Sect. 3. We discuss how the results relate to 
the dynamics of Alfven waves in the solar wind in Sect. 4 
and we summarise our conclusions in Sect. 5. 

2. Mathematical formulation 

2.1. Fundamental properties of Alfven waves 

The equations of ideal isothermal MHD in a homogeneous 
medium can be written as 



| + V.(pv) = 0, 



dt 



V • (vpv) = — Vp + J x B, 



dB 

~dt 



V x (v x B), 



(1) 
(2) 
(3) 
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Fig. 2. Magnetohydrodynamic waves propagating through a homogeneous medium of low j3 and 77 = 0.07 (Run lb), 
a, c, e, g Ap/p — p ~p° versus z/X at t/P — 16.7, 23.3, 33.3, 46.7. Note the changes in the Ap/p-scale between the 
frames, b, d, f, h z_/«a (solid line) and z + /va (dotted line) versus z/X at the same times. 



V - B = 0, 



(4) 



where p is the density, v the velocity, p the pressure, B 
the magnetic field, and J = V x B//io the current den- 
sity. Eqs. 0) and J5J express the conservation of mass 
and momentum, respectively, and Eq. J3J is the induction 
equation. The constraint V • B = is fulfilled by Eq. @ 
if it is imposed as an initial condition. For fully ionised 
hydrogen we can write the equation of state as 



P 



2pk B T 
m H 



(5) 



where fee is Boltzmann's constant, T the temperature and 
mH the mass of the hydrogen atom. 

In a homogeneous medium with a density po and a 
background magnetic field B = B z z a circularly po- 
larised forward propagating Alfven wave is described by 
the transverse magnetic field 



B^ = B±[cos(kz — uit)x + sin(fcz — wt)y], 



(G) 



and the velocity 

v, = --§!==, (7) 

while for a backward propagating Alfven wave the velocity 
is given by 

B, 



(8) 



VM0P0 

In both cases the Alfven wave obeys the dispersion relation 
lo = v^k (9) 
with the Alfven velocity 
B, 



VA 



(10) 



In order to separate forward and backward propagating 
Alfven waves it is useful to introduce the Elsasser variables 



Z+ = 



B 



(11) 
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Table 1. Simulations of Alfven waves in a homogeneous 
medium. The wave is characterised by the two quantities 
77, the amplitude of the imposed Alfven wave in terms of 
the vertical magnetic field, and (3 = 2p p/B 2 the plasma 
beta. For the different Runs we further specify the length 
of the time step At in terms of the period of the wave, P, 
the length of the computational domain, L in terms of the 
wave length, A, and the number of grid points, N. 



Run 


At/P 


N 


L/X 


V 


P 


la 


0.003 3 


3 600 


78 


0.007 


0.042 


lb 


0.003 3 


3 600 


78 


0.07 


0.042 


lc 


0.003 3 


3 600 


78 


0.7 


0.042 


2a 


0.003 3 


13 500 


193 


0.007 


0.96 


2b 


0.003 3 


13 500 


193 


0.07 


0.96 


2c 


0.003 3 


13 500 


193 


0.7 


0.96 


3a 


0.003 3 


13 500 


278 


0.007 


2 


3b 


0.003 3 


13 500 


278 


0.007 


2 


3c 


0.003 3 


13 500 


278 


0.07 


2 


4a 


0.0017 


36 000 


77.2 


0.007 


17 


4b 


0.001 7 


36 000 


77.2 


0.07 


17 


4c 


0.0017 


36 000 


77.2 


0.7 


17 



Galeev & Oraevskii l|1963[) and Sagdeev & Galeev 
(1969) showed that an Alfven wave with a frequency 
w o = va^o and a wave number k Q can decay into a back- 
ward propagating Alfven wave with a frequency w_ and a 
wave number fe_ and a forward propagating acoustic wave 
with a frequency u> and a wave number k that fulfill the 
resonance conditions 



UJ = uj+ I w_ I 
and 

kg = k — k- 



(16) 



(17) 



(e.g. Cramer USD- 

Goldstein l|1978(l and Derby (|1978[l derived the disper- 
sion relation for low frequency waves neglecting dispersive 
effects 



{uj 2 - c 2 k 2 )(w - v A k)[{w + v A k) - Alu 2 



[uj 3 + uj 2 v A k — 3ujuj 2 + ujQV A k] 



(18) 



which describe the forward and backward propagating 
Alfven waves, respectively. For a pure forward propagating 
Alfven wave z- does not exist while z+ has twice the value 
of v± of the original wave. For the circularly polarised 
wave the magnetic pressure B± 2 / (2/io) is the same every- 
where inside the wave. This is the physical reason why an 
infinitely extended circularly polarised Alfven wave is in- 
compressible and an exact solution of the nonlinear MHD 
equations, whereas the linearly polarised wave 

Bj_ = B±[cos(kz — uit)x + cos(fcz — ujt)y], (12) 

is compressible to second order, and therefore not an exact 
solution of the nonlinear MHD equations. 

In the one-dimensional problem that we study, there is 
only one additional wave mode, an acoustic wave obeying 
the dispersion relation 

uj = c s fc, (13) 
where the isothermal speed of sound is 

(14) 

The amplitude of the density oscillation Ap = p — pa is 
related to that of the longitudinal velocity, v z , through 
Ap 

p 

(e.g. Landau & Lifshitz lTW|) . 



(15) 



2.2. Parametric decay 

In the presence of a density fluctuation a circularly po- 
larised Alfven wave decays into forward propagating den- 
sity and magnetic waves in addition to a backward prop- 
agating magnetic wave. These waves are not necessarily 
normal modes of the plasma. 



where 7? = In the limit of /3 <SC 1 Eq. (|T%j) yields a 
maximum growth rate for the instability near k ~ 2fco 
for the density wave. The backward propagating mag- 
netic wave has then w_ ~ — loq and fc_ ~ fco resulting 
in uj- ~ —VAk-. This means that it is a natural mode of 
the system, an Alfven wave, and at the same time the for- 
ward propagating density wave is an acoustic wave. In this 
regime we therefore obtain a parametric instability. When 
(3 increases the daughter wave becomes less resonant with 
the normal mode and the modulational instability loses 
its parametric nature. 

3. Results 

We use the numerical code of Boynton & Torkelsson (1996 ) 
to simulate circularly polarised Alfven waves. This code is 
based on the ETBFCT algorithm fBoris ll1J75|l . a form of 
flux-corrected transport (Boris & Book 1973). The scheme 
is fairly inefficient compared to modern methods; it needs 
about 30 to 40 grid points per wave length of the Alfven 
wave, while there are schemes that require only a third 
of this. However the memory and time requirements for a 
one-dimensional simulation are sufficiently small that this 
is acceptable. 

The waves are driven on the lower boundary of a one- 
dimensional box and propagate through a homogeneous 
medium. In all the runs the grid is sufficiently extended 
that the wave does not hit the upper boundary during 
the course of the simulation. The lower boundary condi- 
tion is formulated in such a way that it only allows an 
outward propagating Alfven wave. This can be a problem 
since backward propagating Alfven waves are generated 
in the simulations. These are mostly of low amplitude 
though, and it is only very late in the simulations that 
high-amplitude Alfven waves reach the lower boundary, 
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Fig. 3. Magnetohydrodynamic waves propagating 
through a homogeneous medium of /3 = 0.96 and 
r) = 0.007 (Run 2a). a Ap/p = versus z/X at 

t/P = 16.7. b fl x /B z (solid line) and v x /v A (dotted 
line) versus z/Xq at the same time, c z + /va (solid line) , 
Z-jvA (dotted line) as versus z/Xq at the same time. 



the amplitude of the Alfven wave. One should note here 
that for a linearly polarised wave the corresponding en- 
hancement shows a strong sinusoidal modulation (1996). 

The density discontinuity at the front of the Alfven 
wave excites a secondary compressional wave which be- 
comes more pronounced for the higher amplitude run lb 
(?/ = 0.07). This wave can be seen as a weak modula- 
tion of the density enhancement (Figs. HJi, Efe), and also 
of the magnetic pressure |B^| 2 /(2/Lto). The density peaks 
coincide with the peaks of the magnetic pressure making 
them a third order analogue of the fluctuations found by 
Boynton & Torkelsson l (T35g| Fig. 2a). 

The density fluctuations serve as the necessary seeds 
of the parametric instability, but due to the low amplitude 
of the fluctuations the instability grows slowly. The evo- 
lution of the backward-propagating Alfven wave that is 
generated by the instability can be followed in Figs. 03, d, 
f and h. It grows in amplitude away from the wave front, 
which enhances the growth rate of the parametric instabil- 
ity upstream. Eventually, the instability is so strong that 
it becomes an efficient source of a forward-propagating 
sound wave (Figs. c, e and g). We get a region at 
z < 15Ao in Fig.[2g and h, in which the waves are strongly 
interacting. In this region there is a strong damping of the 
forward-propagating Alfven wave, and the acoustic wave 
is amplified until it becomes so nonlinear that it steep- 
ens into shocks. On the other hand since the Alfven speed 
is larger than the sound speed, the head of the Alfven 
wave manages to stay ahead of the density fluctuations 
and therefore remains essentially unaffected by the para- 
metric decay. 

In agreement with the theoretical prediction the sound 
wave has a wave number k ~ 2ko. Figures and f show 
that there is a phase shift of ? between the backward- 
and the forward-propagating Alfven waves as we expect if 
the backward-propagating wave is generated through the 
parametric instability. 



and the physical processes that we are interested in occur 
far from the boundary then. 

We study waves of different amplitudes, ry, at different 
values of (3, the ratio of gas pressure to magnetic pressure 
(Tab. [TJ. The magnetic field is strong in runs la-c, that 
is i>a > c s , intermediate in runs 2a-c and 3a-c, that is 
«A ~ c s and weak in runs 4a-c, that is c s > va- 

3.1. Very low j3 

As an example we look at run la (/3 = 0.042) with a 
strong magnetic field and a low-amplitude Alfven wave 
(77 = 0.007). The propagating wave is shown in Fig. QJj. 
The Alfven wave excites a density enhancement (Fig.^Ji), 
whose right edge coincides with the front of the Alfven 
wave, while the left edge propagates with the lower speed 
c s . This density enhancement is a second order effect, and 
its amplitude is consequently proportional to the square of 



3.2. (3 = 0.96 and 2 

We look first at the low amplitude wave 2a. The Alfven 
wave is reflected against a density enhancement at the 
front of the Alfven wave (Fig.[3Ji and c). The density en- 
hancement acts like a wall that is pushed to the right by 
the wave. Fig. |2t shows that the reflected wave, repre- 
sented by 2_, has a phase shift of ir radians relative to 
the mother wave, which is characteristic of reflection at 
a fixed end. To the left of the density enhancement (Fig. 
|3Ji) we see rapid fluctuations in the density (cf. the low /3 
models). 

The right edge of the density enhancement travels at 
a speed ~ c s + v z as expected from an acoustic wave, 
while the left edge of the density enhancement, that is the 
front of the Alfven wave, travels at Vf = 0.23 v a = vx/n. 
As a consequence of this, the measured wave length of 
the Alfven wave is Xq/ti = 0.23 Ao, where Ao = va P- We 
can see also that the wave does not show the ordinary re- 
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Fig. 4. Magnetohydrodynamic waves propagating through a homogeneous medium of j3 = 0.96 and rj — 0.07 (Run 
2b). The plots focus on the Alfven wave front at two times, 16. 6P (a, c) and 66. 6P (b, d). The upper panels a and 
b show Ap/p versus z/Xq, while the lower panels c and d show B x /B z versus z/Xq. 



lation between v x and B x that we expect for an Alfven 
wave (Fig. 1^)). While the transverse velocity component 
maintain the relation v x /va = T) we find for the fluctu- 
ating magnetic field B x /B z = 4.25 t; = rvq. This can be 
easily understood. During half a period the injected flux 
of the x-component of the magnetic field is t]vaB z P/it, 
but because the wave propagates at the speed va/ti, this 
magnetic flux is compressed into a distance vaP/ti, which 
leads to that the magnetic field is amplified by the factor 
n. 

The dynamics of the reflection zone is more easily stud- 
ied in a wave with a higher amplitude, because the extent 
of the reflection zone increases with the amplitude of the 
wave. We show the reflection zone of model 2b at two 
different times in Fig. 0] The density changes smoothly 
along a positive slope between the positions A and B. The 
slope becomes less steep with time since it is spreading out 
over a longer distance. To the left of A the phase speed 
Up = X/P = 0.24 1;^, which is the speed at which the wave 
front is propagating, but between A and B the phase speed 
increases to va- The reason why the wave in run 2b runs 
slightly faster than the wave in run 2a is that it is gaining 
the speed v z of the medium in which it is propagating; a 
similar effect was seen in Fig. 13 of Torkelsson & Boynton 
(1998). 

There is a gradual damping of the Alfven wave to the 
left of A (Fig. EJ- This decay is accompanied by a gradual 



increase of the density. These changes are negligible in the 
low amplitude runs 2a and 3a. 

The higher (3 models 3a and 3b share the characteris- 
tics of models 2a and 2b (with n = 5) , but the damping 
becomes more pronounced. Because of the increase in the 
factor n, the wave in model 3b is slower than the wave 
in run 2b (Fig. |HJ). As shown in Fig. 03 run 2b has lost 
~ 13.3% of its Poynting flux between z = and point A 
at t = 133P. At the same time run 3b has lost ~ 26.5% 
of its Poynting flux. The stronger damping is associated 
with a significant modulation of the Alfven wave on long 
length scales (Fig. EJ 

3.3. = 17 

At a low magnetic pressure (runs 4a-c) the behaviour of 
the wave is almost independent of its amplitude. As an ex- 
ample we take the high amplitude model 4c with r\ = 0.7. 
The speed of the front of the Alfven wave is ~ IAIva 
in models 4a and 4b and ~ 1.43 Va in model 4c. The 
higher wave speed in run 4c is the result of that the Alfven 
wave has generated an outflow with a speed ~ 0.02 va- 
Analogously to models 2 and 3, the amplitude of the mag- 
netic field has changed with a factor 0.7 = to compen- 
sate for the increase in the speed of the wave front (Fig. 
|Eb). The change from sub-Alfvenic waves in runs 2 and 
3 to super- Alfvenic waves in runs 4 is accompanied by a 
change in the character of the reflection. Fig. [7| shows that 
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Fig. 5. Magnetohydrodynamic waves with rj = 0.07 prop- 
agating through a homogeneous medium of (3 — 0.96 (Run 
2b). a Ap/p = versus z/X at t/P = 133. b Poynting 
flux, S Zl versus z/Xo at the same time. So is defined as 
v A B z 2 /hq 
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Fig. 6. Magnetohydrodynamic waves with 77 = 0.07 prop- 
agating through a homogeneous medium of a f3 = 0.96 
(Run 2b) and b f3 = 2 (Run 3b) versus z/X at t/P = 66.7. 



there is no phase delay between the two Elsasser variables 
z± , which is a characteristic of the reflection of a wave at 
a free end. 

Fig. |S| shows the behavior of the wave close to the 
reflection zone. As in models 2 and 3, the slope in density 
becomes less steep with time (compare Fig. [7^ and Fig. 
IHt)- The behaviour of the waves within this zone, in terms 
of the amplitude and the phase speed is similar to that 
described in sect. 3.2 for models 2 and 3 (with n = 1.43). 
The difference is that the damping of the Alfven wave 
outside the reflection zone is too small to be measured. 

4. Discussion 

In constructing models of Alfven wave driven stellar 
winds, it is important to understand the mechanism by 
which the Alfven waves are damped. A linearly polarised 
Alfven wave of high amplitude can steepen and form cur- 
rent sheets even in a homogeneous medium, which leads 
to a quick damping of the wave. While circularly polarised 
Alfven waves cannot steepen in this way, they can be sub- 
ject to a parametric decay into a backward propagating 
magnetic wave and a forward propagating density wave. 
Our simulations show that the compression of the back- 
ground medium that takes place at the wave front of an 
Alfven wave in a low [3 plasma is sufficient to trigger this 
instability. At (3 <; 1 this is not sufficient to trigger the 
parametric decay, but instead we find that the Alfven wave 



is reflected at the wave front, whose speed deviates from 
the Alfven speed. 

One can find a large range of /3s in the solar corona and 
solar wind. Typically (3 <C 1 in the lower solar corona, and 
(3 > 1 in the solar wind fGarv EoTiT]) . This means that dif- 
ferent processes may affect Alfven waves in different parts 
of the solar corona/wind. The parametric decay should be 
at work in the lower solar corona, and in the presence of 
sufficiently large fluctuations it may also be at work in 
the solar wind as indicated by the numerical simulations 
of Del Zanna et al. I2QU1J . Therefore either the paramet- 
ric decay or reflection at the wave front may be the source 
of the inward propagating Alfven waves that have been 
detected in the solar corona (e.g. Bavassano et al. 2000). 

An interesting property of the parametric decay is that 
it is relatively weak at the wave front, but the backward 
propagating Alfven wave that is generated serves as a seed 
for the instability further upstream, where the growth rate 
therefore becomes larger. The area close to the wave front 
is therefore not significantly affected by the parametric 
decay, though the decay becomes strong some distance 
behind the wave front. This gives rise to a self-limitation 
of the length of an Alfven wave packet, and it also suggest 
that, at least in a low-/3 plasma, the Alfven wave is the 
least turbulent right behind the wave front. In a high-/3 
plasma on the other hand one may expect to see a gradual 
change in density right behind the wave front, and further 
upstream only oscillations around a constant density. 
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Fig. 7. Magnetohydrodynamic waves with 77 = 0.7 prop- 
agating through a homogeneous medium of f3 = 17 (Run 
4c). a Ap/p = £z£a versus z/A at t/P = 33. b z + /tM 
(solid line), z^/va (dotted line) and v x /va (dashed line) 
versus z/Ao at the same time. 

5. Conclusions 

In this paper we have presented numerical simulations of 
Alfven waves in homogeneous media with different mag- 
netic field strengths. In a strongly magnetised plasma the 
Alfven wave decays parametrically to a backward prop- 
agating Alfven wave and a forward propagating acoustic 
wave unless the amplitude of the Alfven wave is very low. 
However by its very nature the parametric decay is fairly 
inefficient at the wave front, even when it is the cause of 
a strong damping upstream. 

For (3 ~ 1 or higher we find a different pattern. The 
wave front is propagating at the speed va/ti and the mag- 
netic field fluctuations are amplified by a factor n. n > 1 
for (3 ~ 1 but n < 1 for (3 3> 1. At the wave front the 
Alfven wave is reflected by an extended region with a pos- 
itive density gradient. The spatial extent of this region in- 
creases with time. Inside this region the phase speed of the 
Alfven wave is increasing from va/ti to va- When n > 1 
the reflected Alfven wave is phase-shifted by n, while it 
does not suffer any phase-shift for n < 1. There is also 
a significant gradual damping, <~ 10 — 20% of the Alfven 
wave when n > 1. 

The results presented in this paper are derived from a 
highly simplified model, which allows us to isolate a few 
physical effects and study them in detail. This study can 
guide us in future investigations of the dynamics of Alfven 
waves in more realistic configurations. In forthcoming pa- 
pers we will discuss the effects of stratification and an 
expanding magnetic field. 



Fig. 8. Magnetohydrodynamic waves with rj = 0.7 prop- 
agating through a homogeneous medium of (3 = 17 (Run 
4c). a Ap/p = versus z/\ at t/P = 169. b 



(B x /B z )/r) versus z/X at the same time. 
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